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ABSTRACT 

C I We study the influence of the various parameters of scenarios of large-scale structure 

^ . formation on properties of galaxy clusters, and investigate which cluster properties are 

most sensitive to these parameters. We present a set of large N-body simulations and 
I derive the intrinsic properties of galaxy clusters in these simulations, which represent 

i a volume of 256^ Mpc'^. The cosmological scenarios studied differ in either the 

shape of the power spectrum of initial fluctuations, its normalization, the density 
parameter Qq, or the Hubble parameter Hq. Between each of the simulations, only 
one parameter is set differently, so that we can study the influence of that parameter 
I on the cluster properties. The cluster properties that are studied are the mass, line- 

, of-sight velocity dispersion, peculiar velocity, intrinsic shape, and orientation with 

If) • respect to its surroundings. 

I The present-day r.m.s. mass fluctuation on scales of 8/i~^ Mpc, fjg, which is largely 

. determined by the normalization of the initial power spectrum, has a large impact on 

' the cluster properties. The latter, viz. the cluster number density, mass, line-of-sight 

f-H ^ velocity dispersion and peculiar velocity, are also determined by Jloj though somewhat 

less. Other parameters, such as Hq, the tilt of the initial fluctuation spectrum, and 
the exact shape of this spectrum, are generally less important. 

Unlike the other cluster properties studied, the peculiar velocity is found to depend 
on all parameters of the formation scenario. 
^1 In a companion paper the properties of the model clusters are compared to obser- 

vations to try and discriminate between different cosmological scenarios. Using scaling 
relations between the average properties of the cluster sample and the parameters of 
the formation scenario, one may try and interpolate between the scenarios studied 
here in order to find the parameters of the scenario that is most consistent with the 



o 



' data. 



Key viTords: Galaxies: clustering - cosmology: observations - large-scale structure 
in the Universe. 



1 INTRODUCTION 

Clusters of galaxies are the largest objects in the Universe 
that have recently collapsed, which makes them readily iden- 
tifiable objects whose statistical properties can help to con- 
strain scenarios for large-scale structure formation in the 
Universe. 

N-body simulations of large patches of the Universe 
were used extensively to compare the properties of the sim- 
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ulated clusters with observations. Frenk et al. (1990) cal- 
culated the distribution of velocity dispersions of clusters, 
identified both in 3-D and in 2-D. Jing & Fang (1994) de- 
termined the cumulative distributions of mass, velocity dis- 
persion, temperature of clusters as well as their space density 
in three scenarios, viz. standard CDM, low-fio CDM, and a 
hybrid of CDM and HDM. Jing & Borner (1995) determined 
the velocity dispersion profiles of clusters for seven cosmo- 
logical scenarios. They found that the velocity dispersion 
profiles depend both on Slo and on the cosmological constant 
A. The profiles are steeper in a low-fio scenario than for a 
high value of f2o- This effect is weaker for a larger value of 
A. Mohr et al. (1995) carried out numerical simulations for 
three cosmological scenarios (f2o = 1-0, f2o = 0.2, f2o = 0.2 
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with Ao = A/3Hq — 0.8) with an effective power spectrum 
P{k) oc on cluster scales. They compared the projected 
X-ray shapes of the clusters in these scenarios to the ob- 
served shapes. They favoured the flo ~ 1 scenario, although 
the authors admitted that some discrepancies remain. De 
Theije et al. (1995) used the galaxy distribution of 99 clus- 
ters to determine the projected ellipticities, and compared 
these to the results of N-body simulations of van Kampen 
(1994). They found that in the Q.o = 1 CDM scenario clus- 
ters are too elongated with respect to real observed clusters. 
From a limited number of simulations for an f2o = 0.2 CDM 
model, they concluded that clusters in this scenario are gen- 
erally more nearly spherical than for f2o = 1.0. 

N-body simulations were also used to study cluster 
alignments. Dekel, West & Aarseth (1984) and West, Dekel 
& Oemler (1989) concluded that the orientations of clusters 
of galaxies with respect to their neighbours provide a sen- 
sitive test for the formation of large-scale structure in the 
Universe. 

Many observational properties of clusters were studied 
in the recent past and much numerical work has been done. 
Most of the latter, however, concentrated on just one or two 
cluster properties. Also, in most studies up till now clus- 
ters were identified in the simulations in three dimensions, 
or background galaxies were not removed from the observa- 
tions. 

In order to improve upon this state of affairs, we have 
run a set of cosmological N-body simulations with the fol- 
lowing goals: (1) Evaluate the intrinsic properties of the clus- 
ters, i.e., of the groups of particles defined in 3-D that fulfill 
certain constraints. (2) Investigate the infiuence of the cos- 
mological parameters on the cluster properties and study 
which cluster properties are mostly affected by varying the 
cosmology. (3) Obtain scaling relations between the aver- 
age properties of the cluster sample and the parameters of 
the formation scenario to allow interpolation between the 
studied scenarios. (4) Compare the cluster properties of the 
model clusters with observations. The results of this com- 
parison, together with the results of other studies, e.g., of the 
fiuctuations in the 3K-background measured by COBE, can 
then be used to find the parameters of the scenario that is 
most consistent with all data. For the best-fitting scenario, a 
cluster catalogue will be constructed later on, similar to the 
standard CDM catalogue of van Kampen & Katgert (1997), 
in which each cluster will be simulated individually at higher 
resolution. 

The cosmological parameters that are varied between 
the different scenarios are the shape of the initial power spec- 
trum, its amplitude which gives rise to the present r.m.s. 
mass fluctuation on scales of 8/i~^ Mpc, as, the density in 
units of the critical density, f2o, and the Hubble-parameter, 
Hq. The values of the parameters are chosen such that very 
often one can compare two scenarios that differ by only one 
parametervalue, so that the infiuence of that specific param- 
eter on the cluster properties can be investigated. 

The comparison with observations is done in a compan- 
ion paper (de Theije, van Kampen & Slijkhuis 1997, here- 
after Paper II). In that paper, we will try to mimic as closely 
as possible the observational way of defining a cluster, to 
make a reliable comparison. 

The present paper is organized as follows: in Section ^ 
we describe the N-body simulations, and the identification 



Table 1 . Description of the cosmological scenarios that have been 
examined. The first column gives the acronym, which will be used 
throughout this paper to identify a scenario. SCDM indicates the 
Standard CDM scenario with Oq = 1-0 a^nd ^^ = 0.5. The various 
LCDM scenarios all have Oq < 1.0. hCDM denotes a scenario 
with h = 0.3, and TCDM is a tilted CDM scenario, having a 
slightly tilted initial spectrum. The second column gives the shape 
of the power spectrum. The CDM power spectrum is taken from 
Davis et al. (1985), the HDM one from Bardeen et al. (1986). Col- 
umn 3 denotes the present— day r.m.s. mass fluctuation in spheres 
of 8/i~^ Mpc. For the £7o = 1.0 scenarios, the values are the max- 
imum values we investigate (see text). Column 4 gives the density 
parameter. Column 5 gives the Hubble— constant in units of 100 
km s~^ Mpc~^. Column 6 gives the index of the primordial spec- 
trum, where n = 1.0 denotes the Harrison-Zeldovich spectrum 
and n = 0.8 corresponds to a slightly tilted initial spectrum. 
Column 7 gives the implied present age of the Universe for each 
scenario. 
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1.3 
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0.5 


1.0 
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0.90 


0.2 


0.5 


1.0 


1.7 


LCDMc 


CDM 


0.90 


0.2 


1.0 


1.0 


0.8 


LCDMd 


CDM 


0.60 


0.8 


0.5 


1.0 


1.4 






0.64 


1.0 


0.5 


-2.0 


1.3 


hCDM 


CDM 


1.00 


1.0 


0.3 


1.0 


2.2 


TCDM 


CDM 


0.60 


1.0 


0.5 


0.8 


1.3 


HDM 


HDM 


1.00 


1.0 


1.0 


1.0 


0.7 



of groups in the simulations. In Section |3| we give the num- 
ber densities of clusters in the simulation boxes and compare 
these to the observations. In Section ^ we describe the in- 
trinsic properties of the clusters in the simulations, viz. their 
mass, velocity dispersion, peculiar velocity, shape, and align- 
ment with the surroundings. In Section H we compare the 
properties of clusters for different parameter values, which 
allows us to investigate the influence of a given parameter on 
the cluster properties. Finally, in Section H the main results 
are summarized and discussed. 



2 METHODS 
2.1 Simulations 

The simulations are performed using a P^M code 
(Bertschinger & Gelb 1991) in a box of 256^ h'^ Mpc^ with 
periodic boundary conditions and 128"^ mesh points. 128"^ 
particles are used, each with a mass of 2.22X IO^^^o/i'^Mq. 
Each simulation took about 25 hours of CPU-time on the 
Cray C90 of the SARA Computing Center in Amsterdam. 
A softening length of 0.2 Mpc (i.e., one tenth of a grid 
cell) is used. This value is the same as Frenk et al. (1990) 
adopted, and makes a comparison between both studies pos- 
sible. The accuracy of the simulations is judged from the 
Layzer-Irvine cosmic energy equation (see Efstathiou et al. 
1985). In all simulations, the integration constant of this 
equation changes less than 0.06% of the current potential 
energy at any stage of a run. 

The parameters that define the simulations are given 
in Table hi. All simulations are done using the same ran- 
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dom number seed to set up the initial conditions in order 
to remove cosmic variance from the comparison of the vari- 
ous scenarios. The initial conditions are evaluated at 2; = 9 
by means of the Zeldovich approximation. For the CDM 
simulations we use the Davis et al. (1985) power spectrum. 
For the HDM simulation we use the Bardeen et al. (1986) 
power spectrum. f2o is the present-day value of the density 
parameter, while h describes the present-day value of the 
Hubble-parameter Hq in units of 100 km s~^ Mpc~^. 

In general, erg increases with time, so it is an indication 
of time. For scenarios with Q.q = 1.0, a% evolves linearly 
with the expansion factor (e.g., Padmanabhan 1993) because 
the shape of the power spectrum is constant in time. In 
such a scenario the particle distribution at an epoch before 
the present time is identical to that at the present time 
in a scenario with a correspondingly smaller value of ag. 
Consequently, the values of og, given for the five f2o = 1-0 
scenarios, are the maximum values of erg that one can probe. 
The erg values in Table |l| are roughly in agreement with 
the findings of Eke, Cole & Frenk (1996). These authors 
derived values for erg equal to 0.50 ± 0.04 for Q.o — 1.0 
scenarios, 0.55 ± 0.04 for Qq = 0.8, and 1.03 ± 0.08 for = 
0.2 scenarios, by comparing their N-body simulations with 
the observed number density of clusters as a function of X- 
ray temperature. The LCDMa scenario has a normalization 
similar to the COBE-normalization, erg = 0.46 (Sugiyama 
1995). 

For the scenarios with Qo < 1-0, the shape of the power 
spectrum does change with time and therefore the simple 
rescaling of the time coordinate in the flo = 1.0 scenarios 
is impossible. In these scenarios the particle distribution at 
an epoch before the present time is identical to that at the 
present time in a scenario with a smaller value of as as 
well as larger values of Q and the spectral parameter P. 
The latter was introduced by Efstathiou, Bond & White 
(1992) and measures the spectral shape. For CDM scenarios 
it is equal to O/i. Scenarios with a larger value of P contain 
less power on large scales and slightly more power on small 
scales. For example, for the LCDMa and LCDMb scenario, 
the range in P is 0.10-0.27, for the LCDMc scenario 0.20- 
0.53, and for the LCDMd scenario 0.40-0.47, for expansion 
factors between 1.00 and 0.22. 

n is the index of the primordial spectrum, with the 
canonical value of 1.0 corresponding to the Harrison- 
Zeldovich spectrum. The TCDM scenario has n — 0.8, i.e., 
its spectrum is slightly tilted and has somewhat more power 
on larger scales than standard CDM. 

The LCDMd scenario is included to investigate the dif- 
ference between a scenario with a value of Qq less than but 
close to 1.0 and a scenario with Qq = 1-0. Because of the 
different dynamics between fiat and open scenarios, the fact 
that Qq is smaller than 1.0 may be more important for the 
cluster properties than how much Qo is exactly below 1.0. 
As most studies for Qo < 1-0 so far usually adopt a 'very 
low' value for Qo, i-e., ~ 0.3, this has never been checked 
thoroughly. 

The implied present age of the Universe in each of the 
scenarios is given in column 7 of Table |l]. 

The expansion factor aexp for which the particle posi- 
tions and velocities are stored are almost always aexp ~ 
(0.22, 0.33, 0.44, 0.55, 0.66, 0.79, 0.89, 1.00), except for 
the SCDM and scenarios. For the SCDM scenario we 



have stored aexp = (0.25, 0.33, 0.41, 0.49, 0.59, 0.67, 0.75, 
0.84, 0.90, 1.00) and for the fc"^ scenario aexp = (0.16, 
0.24, 0.31, 0.39, 0.47, 0.56, 0.64, 0.71, 0.86, 1.00). The red- 
shift corresponding to a given expansion factor follows from 
^ = '^exp ~ 1- 

2.2 Defining clusters 

2.2.1 Group finding algorithm 

The algorithm that we used to define groups is the "friends- 
of-friends" algorithm. It is described in, e.g., Davis et al. 
(1985). This algorithm links all particle pairs that are sep- 
arated by less than a fraction p of the mean interparticle 
distance. Each subset of linked particles is then defined as 
a group. The algorithm finds groups that have an overden- 
sity p"'' with respect to the mean background density. Since 
the overdensity Sp/p within the virial radius of a cluster is 
~ 180 (e.g., Padmanabhan 1993), typical values for p should 
be about 0.15-0.20. We adopt p = 0.20. The algorithm has 
the advantage that it produces a unique catalog of groups 
for any p, and that it does not make a priori assumptions 
about the shape of the groups. 

Only groups containing at least 100 particles are in- 
cluded in the cluster list. This lower limit is adopted since 
fewer particles may result in a large increase in shot noise in 
the measurements. In the simulations 100 particles represent 
a mass of 2.22 X IO^/i-^J^oMq. This mass can be compared 
to the mass within the virialized part of clusters, which is 
about 5 X lO"/i~^Af0. 

2.2.2 Cluster definition 

For the study of intrinsic properties of the model clusters, 
all groups with a mass larger than 2.22 x lO^'*/i~^M0 will be 
considered 'clusters' from now on. These clusters form our 
mass-limited cluster catalogue. Each cluster then contains 
at least lOO/ilo particles. This difference in mass resolution 
between Qq ~ 1-0 and flo = 0.2 may influence the analyses. 
In Paper II we will address the problem of mass resolution 
more thoroughly by explicitly undersampling the clusters. 
We conclude there that this undersampling does not have 
any significant infiuence on the cluster properties. 

The observed number density of rich {R > 1) Abell 
clusters is 8.6 x 10"'' /i^ Mpc"^ (e.g., Mazure et al. 1996). 
Thus, one expects about 144 rich clusters in our simulation 
volume. This number of 144 will be used when comparing 
the observed properties of the model clusters to real obser- 
vations (Paper II). As in some of the scenarios the number 
of clusters with a mass above the limit is far less than 144, 
we have chosen to define for each scenario a number-selected 
cluster catalogue as well, consisting of the 144 most massive 
clusters. However, the mass ranges of those catalogues can 
be quite different. For example, in the LCDMb scenario the 
mass of the least massive cluster in the number-selected cat- 
alogues is 1.5 X W^^h-^MQ while for the HDM scenario it 
is 7.7 X W^'^h~^ Mq. The mass of the most massive cluster 
ranges from 5.8 x 10^*/i"^Mq to 2.0 x lO^^/i'^M©. 

For both definitions of the cluster sample, Figure ^ 
shows the percentage of all particles that is in clusters. This 
gives an idea of the importance of such clusters in the simu- 
lations. For the mass-limited samples, two scenarios clearly 
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Figure 1. Percentage of particles in clusters, a: Clusters are 
defined as those groups which have a mass of at least 2.22 X 
10^'^h~^ M0. b: Clusters are defined as the 144 most massive 
groups. 



stand out, viz. the hCDM and HDM scenarios (lefthand 
panel). The reason for this is the large value of erg for these 
scenarios (i.e., erg = 1.00, see Table |^). At an expansion fac- 
tor aexp ~ 0.65, where both scenarios have ag ~ 0.65, the 
fraction of particles in clusters is similar to that in other 
scenarios which also have as ~ 0.65. So, while studying the 
clusters that are defined in this way, one is treating a similar 
fraction of the total number of particles in each simulation 
cube. 

The righthand panel of Figure |l] shows the percentage of 
particles in the 144 most massive groups. Again, the hCDM 
and HDM scenarios stand out, but now the two Qq — 0.2 
scenarios LCDMb and LCDMc have an even larger fraction 
of the particles in clusters, almost three times as much as 
for the other scenarios. If one would rescale these scenarios 
to lower (Tg by investigating the scenario at earlier times, 
the density parameter and the spectral parameter F also 
change. So the fraction of particles involved shows a broader 
range than for the mass-limited samples. However, we use 
the number-selected sample only if the results of the mass- 
hmited sample may be influenced by the limited number of 
clusters in the latter sample. 



2.2.3 Cluster centres 

To determine the cluster centre, two methods are used. The 
first one is the same as in de Theije, Katgert & van Kampen 
(1995). First, one calculates the centre of mass of all particles 
in the cluster. Then one defines an aperture of radius 1.0h~^ 
Mpc around this centre and calculates the new centre of all 
particles within this aperture. This procedure is repeated 
until the mass centre does not shift by more than 0.1h~^ 
Mpc. 

In the second approach to define the cluster centre we 
find the particle which has the largest smoothed density of 
particles around it. This density is calculated by smoothing 
the particle distribution around each particle with a Gaus- 
sian distribution of dispersion Rs, which is equal to half 
the average nearest neighbour distance [3/47r(n)]^''^, where 
(n) is the mean number density of particles (see van Kam- 
pen 1995). This definition resembles the X-ray centre of a 
cluster, because the X-ray emission is proportional to the 
electron density squared of the intracluster gas. 



We have checked that the cluster mass is the same for 
both definitions of cluster centre, and expect this to be true 
for the other cluster properties as well. As the first method 
is much less time consuming, we use that method in what 
follows. 



2.3 Galaxy identification 

Different schemes have been used in the literature to iden- 
tify galaxies in N-body simulations (e.g.. White et al. 1987, 
Nolthenius, Klypin & Primack 1994) . In most of the schemes 
one identifies galaxies as the highest peaks above a thresh- 
old given by some bias factor. Van Kampen (1995) de- 
scribed the formation and evolution of galaxies by replac- 
ing a group of particles that is roughly in virial equilibrium 
by a single soft particle with mass, position, velocity and 
softening corresponding to that group. His constrained ran- 
dom field single cluster simulation cubes were much smaller 
than the present ones and each particle had a mass of 
3.5 X 10^" Mq. As the particles in our simulations have a 
mass of 2.22 X W^^floh'^MQ, already the mass of a galaxy, 
it is impossible to apply this algorithm here. 

Instead, we assume the dark matter to be equally dis- 
tributed as the luminous matter, i.e., galaxies and gas. Buote 
& Canizares (1996) concluded that, for their sample of 5 
clusters, the shapes of the dark matter distribution, the mass 
distribution and galaxy isopleths are all consistent with each 
other. In addition, for the catalogue of model clusters of 
van Kampen (1995), the distributions of ellipticities for the 
galaxy particles on the one hand, and for all dark matter 
particles on the other hand are statistically the same. 

From a sample of 41 clusters with measured veloc- 
ity dispersions and X-ray temperatures, Lubin & Bahcall 
(1993) concluded that the velocity bias in clusters, = 
"^gal/^^DM ~ ^■^'^ 0.04. So the velocity dispersion of the 
galaxies is almost identical to that of the dark matter. Also, 
van Kampen (1995) found no evidence for velocity bias in his 
catalogue of cluster models which contain a recipe for galaxy 
formation. In addition, several authors concluded that the 
spatial distribution of the galaxies, gas and the total mass 
are all very similar, with possibly a somewhat larger cen- 
tral concentration for the dark matter (e.g., Henry, Briel 
& Nulsen 1994, Smail et al. 1995, Tyson & Fischer 1995, 
Squires et al. 1996). 



3 CLUSTER NUMBER DENSITIES 

An important property of the cluster distribution as a whole 
is the number density of rich Abell clusters in the simula- 
tions. In the present paper we define the cluster sample to 
be all groups in the simulation that have a mass of at least 
2.22 X lO^'*/i"^M0 (see Section ^.2.2| ), i.e., the mass-limited 
catalogue. 

In Table |^ we give the number of clusters in each sim- 
ulation, i.e., the number of groups which have a mass of 
at least 2.22 x 10^'^h~^ Mq within a clustercentric radius of 
l.Oh"^ Mpc at various times. For the SCDM scenario, these 
values are similar to the predicted number of clusters above 
this minimum mass as derived from the Press-Schechter for- 
malism (Eke et al. 1996). The numbers in parentheses in Ta- 
ble 1^ give the mean number of particles per cluster within 
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Table 2. Number of clusters with a mass within r = l.Oh ^ Mpc of at least 2.22 X 10^** /i ^Mq, corresponding to lOO/Oo particles. The 
numbers in parentheses give the mean number of particles per cluster, averaged over all clusters. 



scenario 



expansion factor 





0.22 


0.33 


0.44 


0.55 


0.66 


0.79 


0.89 


1.00 






LCDMa 






























LCDMb 


1 


5 (319) 


17 (360) 


30 (396) 


29 (501) 


39 (567) 


44 (637) 


46 (679) 






LCDMc 


1 


5 (316) 


17 (341) 


28 (396) 


32 (461) 


42 (563) 


41 (644) 


41 (662) 






LCDMd 











6 (127) 


18 (138) 


51 (146) 


100 (160) 


181 (164) 






hCDM 





2 (108) 


30 (130) 


181 (134) 


424 (143) 


749 (159) 


962 (170) 


1166 (176) 






TCDM 











3 (108) 


18 (131) 


82 (126) 


180 (132) 


338 (136) 






HDM 





2 (121) 


38 (129) 


205 (139) 


455 (151) 


747 (163) 


928 (175) 


1062 (187) 








scenario 










expansion factor 










0.25 


0.33 


0.41 


0.49 


0.59 


0.67 


0.75 


0.84 


0.90 


1.00 


SCDM 











1 (106) 


12 (128) 


34 (122) 


77 (124) 


159 (127) 


232 (134) 


406 (136) 




scenario 










expansion factor 










0.16 


0.24 


0.31 


0.39 


0.47 


0.56 


0.64 


0.71 


0.86 


1.00 



















1 (100) 


8 (131) 


23 (125) 


72 (147) 


189 (149) 



l.Ohr^ Mpc, averaged over all clusters. This mean number 
of particles increases towards later times for all cosmological 
scenarios, illustrating the growth of the clusters. 

In the LCDMa scenario, no clusters with M > 2.22 x 
are formed for the assumed value ag = 0.46 in 
the simulation. Therefore, we will not consider this scenario 
any further. 

If one evaluates the total mass of the clusters, using the 
entire group as defined by the friends-of-friends algorithm, 
the results are qualitatively the same as when one uses the 
mass within 1.0h~^ Mpc. 



4 INTRINSIC PROPERTIES OF CLUSTERS 

The intrinsic properties of the model clusters are evaluated 
using both the mass-limited and the number-selected clus- 
ter catalogues. In the following analyses, all cluster particles 
that are within a clustercentric radius of r = l.Oh^^ Mpc 
are considered, unless stated otherwise. 



4.1 Masses 

The first, most basic, property of the clusters that we will 
investigate is their mass M. Figure ^ shows the cumula- 
tive distribution p(> Af ) of particle mass for all clusters in 
the mass-limited sample. For the Qo ~ 1-0 scenarios, dif- 
ferent values of as are plotted. In these cases, as is highest 
for the upper curves and decreases downward. In the up- 
per righthand panel, the low-no scenarios are shown. For 
the r^o = 1.0 scenarios, the cluster mass increases with as, 
or, equivalently, with time. If one were to extrapolate this 
to as = 0.9, the value adopted for the LCDMb and LCDMc 
scenarios, the clusters in the f2o = 1-0 scenarios have a much 
larger mass than those in these two low-f2o scenarios for 
the same value of as- Even in the LCDMd scenario with 



Slo = 0.8, the cluster mass is significantly smaller than in 
the SCDM scenario with as ~ 0.61. 

A comparison of the distributions for the LCDMb and 
LCDMc scenarios is a check of the calculations, as these 
scenarios are exactly the same apart from the value of the 
Hubble parameter h. As the masses here are expressed in 
units of 10^'^h~^ Mq, there should be no difference at all 
between these two scenarios for their cluster masses. Indeed, 
we find very good agreement. 

Among the different scenarios, the particle mass of clus- 
ters may be quite different. The low-ilo scenarios have 
a most massive cluster of about 6 x 10^*/i~^AfQ. The 
SCDM, and TCDM scenarios contain masses up to 
8 X 10"/i"^A/Q. The hCDM and HDM scenarios con- 
tain clusters of even larger mass, 18 x 10^'^ Mq. This is 
due to the high value of as, viz. erg = 1.00, in these latter 
two scenarios. 

For all flo = 1-0 scenarios the slope of p(> M) is clearly 
flatter at later times, showing that the relative number of 
high-mass clusters increases at later times. The change in 
the median value of the cluster particle mass, however, is 
not very large. This is clear from Table ^, which shows the 
median value of the cluster particle mass for different values 
of aexp. Only for the hCDM and HDM scenarios is there 
an increase in the median value towards later times. For all 
other scenarios the increase is not signiflcant or not mono- 
tonic. 

The median value of the cluster particle mass in the 
low-f2o scenarios changes less rapidly with time. To illus- 
trate the change of cluster mass in these open scenarios, we 
show in Figure ^ the change of the particle mass distribution 
of the low-f2o scenarios LCDMc and LCDMd with time. For 
the LCDMc scenario, the cluster population does not evolve 
strongly for redshifts z < 0.50. I.e., the total number of 
clusters is more or less constant and only a few high-mass 
clusters increase their mass even further. This result is quan- 
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Figure 2. Cumulative distribution of cluster particle masses for 
all clusters having a particle mass of at least 2.22 x lO^^h-^MQ. a: 
SCDM with (Tg = 0.46 (solid line) and erg = 0.61 (dotted line), b: 
LCDMb (solid line), LCDMc (dotted line) and LCDMd (dashed 
line), c: with trg = 0.46 (solid line) and eg = 0.64 (dotted 
line), d: hCDM with eg = 0.44 (solid line), erg = 0.65 (dotted 
line) and erg = 1.00 (dashed line), e: TCDM with erg = 0.47 
(solid line) and erg = 0.60 (dotted line), f: HDM with o-g = 0.44 
(solid line), erg = 0.66 (dotted line) and trg = 1.00 (dashed line). 
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Figure 3. Time evolution of the cumulative distribution of 
cluster particle mass. Only clusters with a mass of at least 
2.22x 10^'*h~^MQ are considered, a: LCDMc scenario. The curves 
correspond to redshifts z = 1.3 (solid line), 0.53 (dotted line), 0.27 
(dashed line), and 0.0 (dot-dashed line), b: LCDMd scenario. The 
curves correspond to redshifts z = 0.53 (solid line), 0.27 (dotted 
line), 0.12 (dashed line), and 0.0 (dot— dashed line). 



Table 3 . Evolution parameters of the mass of all clusters having a 
mass of at least 2.22x10^'^ Mq . The values in this table denote 
the median values of the cluster particle mass M (in units of 
10^'^h~^ Mq). The information is listed only for expansion factors 
at which there are at least 10 clusters. 

median value of cluster particle mass (in W^'^h^^ Mq) 



scenario 



expansion factor 





0.44 


0.55 


0.66 


0.79 


0.89 


1.00 


LCDMb 






2.9 


2.9 


3.0 


2.9 


LCDMc 








2.8 


3.1 


2.8 


LCDMd 






2.6 


2.5 


2.6 


2.7 


hCDM 


2.9 


2.7 


2.8 


3.1 


3.2 


3.3 


TCDM 






2.7 


2.6 


2.7 


2.7 


HDM 


2.7 


2.8 


3.0 


3.2 


3.3 


3.5 





scenario 






expansion factor 






0.59 


0.67 


0.75 0.84 0.90 


1.00 


SCDM 


2.8 


2.7 


2.6 2.6 2.8 


2.8 




scenario 






expansion factor 






0.47 


0.56 


0.64 0.71 0.86 


1.00 








2.7 3.0 


2.8 



titatively consistent with that of, e.g.. Eke et al. (1996). In 
the LCDMd scenario, the situation is quite different. Here 
the number of clusters increases up to the present time and 
the maximum mass increases somewhat faster than in the 
LCDMc scenario. So there is a clear difference in the time- 
evolution of cluster mass with f2o, even between scenarios 
with 0,0 < 1.0. Alternatively, the different curves in Figure 
^ can also be viewed as corresponding to scenarios with a 
different value of F (see Section 2.1). In that case, F is low- 
est for the upper curves. Scenarios with a larger value of F 
then seem to have less clusters above a certain mass, but 
this result is solely due to the different values of erg that 
correspond to the different curves. 

If one uses the average mass estimator (Heisler, 
Tremaine & Bahcall 1985) 

^'^ '^{v^i-VzjfRij, (1) 



Mi 



GN{N- 1) 



where G is the gravitational constant. A'' the total number 
of particles, Vzi is the velocity of galaxy i w.r.t. the clus- 
ter center, and Rij is the distance between galaxies i and 
j, then the cumulative distributions of cluster masses look 
very similar to the distributions of cluster particle mass. The 
average mass estimator, however, overpredicts the cluster 
particle mass by about 40% in all scenarios. The projected 
mass estimator overpredicts the cluster particle mass even 
more than the average mass estimator. These results are 
quantitatively similar to those of van Kampen (1995), who 
used N-body simulations of single clusters in an fig = 1.0 
CDM Universe to test the mass estimators. He found that 
the average mass estimator overpredicts the cluster mass by 
about 40%, if evaluated within about 1.0h~^ Mpc, and that 
the projected mass estimator overpredicts the cluster mass 
by about 60%. 
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Figure 4. Cumulative distribution of line— of— sight velocity dis- 
persion for clusters in the different cosmological scenarios. Only 
clusters with a mass of at least 2.22 X 10^'^h~^ Mq are consid- 
ered. The different panels and curves have the same meaning as 
in Figure 2. 



To investigate how the mass increases with cluster- 
centric radius, we determine the mass of the clusters within 
an Abell radius, 1.5/i~^ Mpc. Using the particle mass for 
this, the mass within the Abell radius is about 25-35% larger 
than that within l.Qh~^ Mpc. This is more or less what is 
expected for a spherical cluster whose density profile follows 
a modified Hubble law p(r) = p(0)/[l + {r/rcff^^ with a 
core radius of rc = 0.25/i~^ Mpc (e.g., Sarazin 1986). 

If one uses the number-selected cluster samples, the 
cumulative distributions of cluster mass would be very sim- 
ilar to the curves in Figure ^ of the mass-limited sam- 
ples. The distributions of cluster average mass change in 
the same way as those over cluster particle mass, but the 
transition between p(> M) < 8.6 x 10"^/i^ Mpc"^ and 
p(> M) > 8.6 X IQ-'^/i^ Mpc"^ is now rather smooth be- 
cause of the scatter in the relation between a cluster's par- 
ticle mass and average mass. 



4.2 Line— of— sight velocity dispersions 

The next property of clusters that we consider is their line- 
of-sight velocity dispersion fj^g, which describes the dy- 
namical state of the cluster, as it is influenced by merging 
events, substructure and the shape of the galaxy orbits. 



Table 4. Evolution parameters of the line-of-sight velocity dis- 
persion of all clusters having a mass of at least 2.22xl0^'^h-'^ Mq. 
The values in this table denote the median values of the line-of- 
sight velocity dispersion crj^g (in km s~^). The information is 
listed only for expansion factors at which there at least 10 clus- 
ters. 

median value of cluster line— of— sight velocity dispersion 



scenario 






expansion factor 




0.44 


0.55 


0.66 


0.79 


0.89 


1.00 


LCDMb 






652 


665 


672 


683 


LCDMc 








654 


657 


635 


LCDMd 






625 


661 


646 


619 


hCDM 


670 


638 


642 


661 


679 


690 


TCDM 






634 


654 


628 


638 


HDM 


663 


649 


669 


686 


702 


719 


scenario 






expansion factor 




0.59 


0.67 


0.75 


0.84 


0.90 


1.00 


SCDM 


568 


650 


645 


632 


638 


637 



expansion factor 



0.47 0.56 0.64 0.71 0.86 



1.00 



637 679 



665 



Crone & Geller (1995) studied the evolution of the clus- 
ter velocity dispersion using large-scale N-body simulations. 
They concluded that crj^g can be significantly influenced by 
merger activity and therefore does not simply refiect the 
cluster mass (they found a scatter of about 5-10% in crj^g 
for clusters of the same mass). Furthermore, they detected 
some change in aj^g with time. The slope of the cumulative 
distribution of the number density of clusters with a velocity 
dispersion larger than crj^g, p(> crj^g), flattens with time, so 
that at later times more clusters are found with large cj^g. 
The change with time is most evident for scenarios with a 
high value of Qq- 

The velocity dispersion for the clusters in our N-body 
simulations is calculated using the robust biweight estimator 
of Beers, Flynn & Gebhardt (1990), yet the velocity disper- 
sions are essentially the same with the ordinary definition. 
The cumulative distributions of line-of-sight velocity dis- 
persion found in our simulations are plotted in Figure ^. As 
before, these distributions include all clusters which have a 
mass of at least 2.22 x 10^* Mq. The different panels and 
curves have the same meaning as in Figure |^, i.e., for the 
rio = 1-0 scenarios the higher curves correspond to larger 
values of erg. 

For the Qq = 1.0 scenarios, our results are consistent 
with those of Jing & Fang (1994) and Crone & Geller (1995). 
That is, the curves for lower erg, or, alternatively, higher z 
(see Section 2.1), are steeper than the z = curves (the 



upper ones). Clearly, the evolution of velocity dispersions is 
not self-similar: the curve shifts towards higher fj^g, but at 
the same time it fiattens. At the present epoch, clusters with 
a high (TjQg are relatively more abundant than clusters with 
a low value of crj^g. These changes in time are analogous to 
that for the mass of the clusters (see the previous section). 
Table ^ lists the median values of o-j^g for the cluster sam- 
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Figure 5. Time evolution of the cumulative distribution of clus- 
ter velocity dispersions. Only clusters with a mass of at least 
2.22x 10^*/i~^Mq are considered, a: LCDMc scenario. The curves 
correspond to redshifts z = 1.3 (solid line), 0.53 (dotted line), 0.27 
(dashed line), and 0.0 (dot-dashed line), b: LCDMd scenario. The 
curves correspond to redshifts z = 0.83 (solid line), 0.53 (dotted 
line), 0.27 (dashed line), and 0.0 (dot— dashed line). 



Figure 6. Scatter plot of the cluster line— of— sight velocity dis- 
persion cr^Qg versus its mass for clusters in the SCDM scenario. 
The solid line shows the relation a oc \/ A/ave . The horizontal and 
vertical dotted lines indicate the values of crj^g and M, respec- 
tively, of the 144-th most massive group, in terms of the particle 
mass, a: crj^g versus cluster particle mass, b: crj^g versus cluster 
average mass estimator Mave. 



pie at all time steps and for all scenarios for which there 
is a significant number of clusters present that satisfy the 
minimum mass criterion. The scenarios with a low erg have 
relatively smaller values of crj^g. This is because less matter 
has collapsed on cluster scales. On the other hand, clusters 
that are not in the sample at early times, because their mass 
is too low, can enter the sample at later times and populate 
the low-cTjQg part of the distribution. The combination of 
both facts results in an apparent upward shift of the whole 
distribution. Only the hCDM and HDM scenarios show a 
significant increase in the median value of crj^g with time. 
For both scenarios we already concluded that the median 
values of th e cluster particle mass increase with time (see 
Section ' 



As an example of the evolution of the velocity disper- 
sions in low-f2o scenarios, we show p(> fj^g) for the LCDMc 
and LCDMd scenarios in Figure]^. For the LCDMc scenario, 
p(> crjQg) hardly changes for z < 0.50. The number of high- 
(TjQg clusters remains constant, as does the maximum value 
of o-jQg. This is consistent with results of Jing & Fang (1994) 
who found no evolution in p(> crj^g) for z < 0.5 for their 
sample of clusters with mass larger than 1.7 x 10^'^ M0. Note 
that these authors used Qq = 0.3. In addition, o-j^g does not 
exceed 1000 km s~^, in contrast to most of the f^o — 1.0 sce- 
narios. For the LCDMd scenario the situation is quite dif- 
ferent. The number of clusters rises until the present time. 
The median value of crj^g, however, is fairly constant. For 
the different values of flo, p(> fj^g) ch ange s with time in a 
similar manner as p(> M) (see Section 4.1). 

For the hCDM and HDM scenarios the fractional 
change in mass is, in general, about twice as large as for 
the velocity dispersion (Tables ^ and ^) , consistent with the 
virial theorem estimates. This evolutionary difference be- 
tween mass and velocity dispersion is also consistent with 
the findings of Crone & Geller (1995), who found that the 
velocity dispersion evolves less rapidly than the mass of a 
cluster. They attributed this to two processes: first, for a 
specific mass range, the velocity dispersion decreases with 
time due to relaxation. Secondly, mergers will make the dis- 
tribution of (TjQg more irregular than that over mass. That 



is, the scatter in the aj^g — M-relation due to mergers and 
accretions introduces random fiuctuations in p(> crj^g). 

To check how large the scatter between (jj^g and M is, 
we show in the lefthand panel of Figure ^ the line-of-sight 
velocity dispersion versus particle mass for all clusters in the 
SCDM scenario for as — 0.61. The solid line shows, for com- 
parison, the best-fitting relation aj^g oc -^/M, as expected 
for systems in virial equilibrium. The r.m.s. scatter about 
this relation is about 14%, somewhat larger than the 5-10% 
that Crone & Geller (1995) found. These results are sim- 
ilar for all scenarios. For comparison, the righthand panel 
of Figure ^ shows the same scatter plot but using the av- 
erage mass estimator. The scatter about the linear relation 
o'lQg oc \/Mave is smaller than if using the cluster particle 
mass, namely 7%. So although the average mass estimator 
overestimates the cluster particle mass by about 40% (see 
Section it correlates better with the cluster line-of- 

sight velocity dispersion than does the cluster particle mass. 
This is because both the velocity dispersion and the average 
mass estimator explicitly use the particle velocities whereas 
the particle mass does not contain this information. 

The horizontal and vertical dotted lines in Figure ^ in- 
dicate the values of aj^g and M, respectively, of the 144-th 
most massive cluster, in terms of the particle mass. These 
values will be called (J144 and M144 from here on. This com- 
parison gives an idea of how complete, in terms of mass, a 
sample of clusters will be that is selected on the basis of their 
velocity dispersion. This is the case, e.g., for the ENACS- 
survey (Katgert et al. 1996), that is claimed to be complete 
for <J\Qg > 800 km s~^. Remember that we expect 144 clus- 
ters in a simulation box of (256 Mpc)'^ on the basis 
of the cluster number density of the ENACS-survey. When 
one uses the cluster particle mass, the fraction of clusters 
with CTjQg > (J144 that also have M > M144 ranges from 64% 
(SCDM) to 78% (for the scenario). These numbers are 
comparable to the 66% of clusters of richness class R > 1 
in the Abell catalogue that also has the required intrinsic 
richness (van Haarlem, Frenk & White 1997). 

To check if the velocity dispersion changes with cluster- 
centric radius, we calculate the value of o-j^g within a clus- 
tercentric radius r = 1.5h~^ Mpc. The velocity dispersions 
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Figure 7. Distributions of 1-Dl peculiar velocities of clusters 
in the different cosmological scenarios. The different panels and 
curves have the same meaning as in Figure 2. The solid lines show 
the best— fitting Gaussian distributions. The parameters of these 
are listed in Table 4. 



within the canonical Abell radius of 1.5/i~^ Mpc are, on av- 
erage, 5% smaller than those within l.Qh~^ Mpc. The exact 
numbers range from 2% to 7% and depend on the scenario 
and on the value of erg. 

For the number-selected cluster samples, the curves of 
P(> '-'^los) would be very similar to those of the mass-limited 
samples. For those scenarios that have at least 144 clusters in 
the mass-limited samples, the curve for the number-selected 
sample is equal to that of the mass-limited sample for p(> 
f^los) ^ 8-6 X lO"*^/!^ Mpc"^, and flattens off smoothly to 
this constant value for lower fj^g. The transition is smooth 
because of the scatter in the o-j^g — M-relation (see Figure |^. 
For scenarios that have less than 144 clusters in the mass- 
hmited sample, the curve for the number-selected sample 
has the same slope for high o-j^g and extends to lower crj^g 
until it flattens off towards p(> aj^g) = 8.6x 10"''/i^ Mpc"^. 

4.3 Peculiar velocities 

The next property of galaxy clusters that we consider is the 
distribution of cluster peculiar velocities, i.e., the velocity 
of a cluster with respect to the Hubble-flow. This peculiar 
velocity is calculated using the robust biweight estimator of 
Beers et al. (1990). 



Table 5. Parameter values of the Gaussian fits to the distribu- 
tions of cluster peculiar velocities. The first two columns specify 
the scenario. The third and fourth columns describe the best- 
fitting Gaussian distribution, /(fpec 1d) ^ '^-''P(~(''pec ID ~ 
'^poc ID o)^/^''^pec ID^' '^^^ fifth column gives the dispersion of 
the Maxwellian fits to the distribution of 3-D cluster peculiar 
velocities. Column 6 gives the predictions of BBKS using linear 
theory. See text for more details. 



scenario crs fpec,lD,0 '^pec,lD '^pec,3D '^pec,lD,lin 







(km s~^) 


(km s"-*^) 


(km s~^) 


(kms ) 


C^ TV T 

bCDM 


0.46 


-16.7 


159 


201 


125 


O T~\ TV T 

bCDM 


0.61 


-4.4 


251 


267 


166 


LCDMb 


0.90 


12.0 


182 


216 


232 


LCDMc 


0.90 


29.4 


160 


181 


156 


LCDMd 


0.60 


-7.3 


223 


247 


162 




0.46 


22.9 


308 


309 




fc-2 


0.64 


-4.6 


447 


508 




hCDM 


0.44 


-39.3 


225 


242 


158 


hCDM 


0.65 


8.9 


292 


318 


234 


hCDM 


1.00 


-4.6 


419 


438 


360 


TCDM 


0.47 


-14.0 


190 


230 


145 


TCDM 


0.60 


-4.2 


262 


280 


186 


HDM 


0.44 


-11.1 


157 


184 


97 


HDM 


0.66 


-18.0 


238 


254 


145 


HDM 


1.00 


-8.0 


333 


341 


221 



Figure [7| shows the differential distributions of 1-D pe- 
culiar velocity I'pgg iy) clusters in the mass-limited 
sample in all scenarios. Also shown are the best-fitting 
Gaussian distributions. Such Gaussian distributions are ex- 
pected for Gaussian random fields (eq. 4.23 of Bardeen et 
al. 1986, hereafter BBKS). In general, the Gaussian distri- 
butions provide good fits. The fitting parameters are given 
in columns 3 and 4 of Table ^ The mean value for the best- 
fitting Gaussian distributions always differ from zero by less 
than 40 km s~^. This is expected because there is neither 
any preferred direction in the simulation box nor any sys- 
tematic fiows on the scale of the simulation box. The disper- 
sion of the fitted Gaussian distributions increases with erg, 
or equivalently cosmic time, for the f2o = 1-0 scenarios. 

For a specific value of erg, the dispersion of the fitted 
Gaussian distribution is smaller for the f2o < 1.0 scenarios 
than for the Ho = 1.0 scenarios, consistent with earlier find- 
ings of, e.g., Bahcall, Gramann & Cen (1994). Again, this 
is because in the open scenarios less matter has collapsed 
on large scales and the total gravitational force acting on a 
cluster is smaller. 

In order to check if the cluster peculiar velocities are 
distributed isotropically, we investigate the distribution of 
3-D peculiar velocities, «pgg 3)3. If the cluster peculiar ve- 
locities are distributed isotropically, the 3-D peculiar ve- 
locity vectors should have a random orientation. We check 
this by determining the distributions of spherical angles cj) 
and cos (6) of the velocity vector. These distributions are 
indeed consistent with uniform distributions for all scenar- 
ios. Column 5 of Table ^ lists the dispersions cr ^ gj-j of 



© 0000 RAS, MNRAS 000, 



10 p. A.M. de Theije, E. van Kampen & R.G. Slijkhuis 



0.50 
0.40 
0.30 
0.20 

0.10 

0.00 

0.40 
0.60 
0.20 

o.-;o 

0.00 



SCDM 

n„ = 1.0 

h = 0.5 



LCDM 

n„ = 0.2,0.8 
h = 0.5,1.0 



n„ = 1.0 

h = 0.5 



I | -r1 
hCDM 

n„ = 1.0 

h = 0.3 



0.40 


. n„ = 1.0 




h = 0.5 


0.30 


. r = 0.8 


0.20 




0.10 




0.00 






0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 1.0 

El El 

Figure 8. Distribution of intrinsic cluster ellipticities ei for dif- 
ferent cosmological scenarios, ei = 0.0 means that the clus- 
ter is spherical, ei = 1.0 means that the particle distribution 
of the cluster is a disk. Only clusters with a mass of at least 
2.22 X 10^*h~^ Mq are taken into account. The different panels 
and curves have the same meaning as in Figure 2. 



the Maxwellian fits /(«pec,3D) °^ "pec,3D exp(-(«pec,3D - 
^pec 3D O^^/^^^pec 3D-^ distributions of 3-D cluster 

peculiar velocity. Column 6 gives the expected 1-D disper- 
sions using linear theory (BBKS) and peaks on a scale of 
4h~^ Mpc. For the scenario, these values cannot be 
obtained because the integral to calculate CTpecSD diverges 
(see e.g. Padmanabhan 1993). The values of cTpg^ for 
the model clusters are almost always larger than those of 
the peaks in linear theory. This indicates that clusters have 
evolved into the non-linear regime in most of the scenarios. 



4.4 Shapes 

De Theije et al. (1995) used the results of N-body simu- 
lations of van Kampen (1994) to study how the shapes of 
clusters may depend on cosmological scenario. They found 
that in a low-f7o CDM Universe, clusters are, on aver- 
age, more nearly spherical than the same clusters in an 
Q.0 = 1.0 Universe. This conclusion was also reached by 
Mohr et al. (1995) who determined cluster morphologies 
from SPH-simulations. These calculations showed that the 
X-ray shapes of clusters are less flattened spherical for low 
values of S7o than in an Einstein-de Sitter Universe. Wil- 



son, Cole & Frenk (1996) also concluded that clusters in an 
Slo = 0.2 Universe are more nearly spherical and centrally 
concentrated than clusters in an Slo = 1-0 Universe. On the 
basis of this result they constructed a new lensing statistic 
that is very sensitive to the value of Q.o and almost indepen- 
dent of the value of the cosmological constant A. 



4.4.1 Distribution of intrinsic shape 

To investigate which cosmological parameters influence the 
shape of a cluster, we determine the shape of the clusters 
for all scenarios. The present simulations are unconstrained, 
unlike the simulations of van Kampen & Katgert (1997), and 
we cannot compare individual clusters in different scenarios. 
Instead, we have to compare the entire distribution of cluster 
shapes. We describe the shape by two ellipticities which are 
obtained using the tensor of inertia, which was shown to 
yield reliable cluster ellipticities (de Theije et al. 1995). The 
tensor of inertia is defined as 



k 



= 1,3), 



(2) 



where the sum is over all particles in the cluster, Xik {i = 
1, 3) are the coordinates of particle k with respect to the clus- 
ter centre and rk is the distance of that particle to the cluster 
centre. The intrinsic cluster ellipticities are then given by 

b 



ei = 1- -, 



£2 



1- -, 



(3) 



where a > b > c are the eigenvalues of the tensor of inertia 
/ (Eq. The triaxiality parameter T was introduced by 
Franx, lUingworth & de Zeeuw (1990): 



T = 



£2(2-62) 
ei(2-ei) 



(4) 



A value of r = 1.0 indicates that a cluster is prolate while a 
value of r = 0.0 represents an oblate cluster. Values between 
0.0 and 1.0 describe triaxial clusters for which a, b and c all 
have different values. 

Figure ^ shows the distributions of ei. There are no 
obvious differences among the different scenarios. For all 
scenarios there is a very small number of spherical clusters 
(ei = 0) and the largest cluster ellipticity is about 0.8. How- 
ever, the Kolmogorov-Smirnov (KS from here on) test (e.g.. 
Press et al. 1989), shows some differences between different 
scenarios. The scenario with as = 0.64 and the HDM 
scenario are significantly different from the SCDM, LCDMb, 
LCDMc, LCDMd, hCDM and TCDM scenarios. Clusters in 
the fc"^ and HDM scenarios have somewhat smaller elliptic- 
ities than those in the other scenarios. The distributions for 
different values of are always consistent with each other 
for all f2o = 1-0 scenarios (the KS-confidence levels for these 
scenarios are always larger than 0.23). 

The ellipticities within the Abell radius (r = 1.5h~^ 
Mpc) are somewhat larger than those within r — 1.0h~^ 
Mpc. The difference is about Aei ~ 0.05 and occurs in all 
scenarios. Qualitatively, such an effect is expected for elon- 
gated clusters: the aperture bias, as a result of which the 
cluster seems more spherical than it actually is (e.g., de 
Theije et al. 1995), is less important for a large aperture 
radius. The projected position ang le is found to change very 
little with radius (see Section [l.5|). 



© 0000 RAS, MNRAS 000, 



Intrinsic Properties of Galaxy Clusters in Different Cosmological Scenarios 11 




0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 1.0 
T T 
Figure 9. Distribution of triaxiality parameter T for clusters in 
the different scenarios. Only the clusters that have a mass of at 
least 2.22 X lQ^^h~^ Mq are taken into account. A value T = 1.0 
means that the cluster is prolate, T = 0.0 indicates that the 
cluster is oblate. The different panels and curves have the same 
meaning as in Figure 2. 



Figure |9| shows the distribution of the triaxiality param- 
eter T for all scenarios. The different curves have the same 
meaning as in Figure ^. No obvious differences are found 
in the distribution of T between different scenarios. In all 
scenarios a large number of clusters have T — 0.6 — 0.9, in- 
dicating that most clusters are nearly prolate. Perfect oblate 
clusters {T = 0.0) are absent in most scenarios. The low- 
no scenarios contain clusters of slightly larger values of T. 
The 'strange' distribution of T for the scenario with 
(Tg = 0.40 is probably the result of limited statistics: the 
number of clusters is only 23 (see Table 



4-4-2 Evolutionary changes in ei, £2 and T 

The distributions of ei for different values of erg in the fio = 
1.0 scenarios in Figure ^ are not significantly different. This 
does not necessarily mean that individual clusters have a 
constant ei and €2- To investigate this, we determine ei and 
t2 for the clusters in the hCDM scenario at different times, 
as an example for the f2o = 1-0 scenarios. In Figures |lo| a and 
b, the values of ei and £2 at 2 = 1.3 are compared with those 
at z = 0.53. For this, only clusters that fulfil the minimum 
mass limit at both redshifts are used. We identify the same 
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e, (z = 0.53) £2 (z = 0.53) 

Figure 10. Evolution of the shape of clusters in the hCDM 
scenario. Only clusters are used that have a mass of at least 
2.22 X IQi^/i-lM© at both redshifts. a: ei at 2 = 1.3 versus 
z = 0.53. b: £2 at 2 = 1.3 versus z = 0.53. c: ei at 2 = 0.53 versus 
2 = 0. c: £2 at 2 = 0.53 versus z = 0. 



clusters at both redshifts by requiring that at least half of 
the particles of the cluster at one redshift is member of the 
cluster at the other redshift, and vice versa. This uniquely 
links clusters at different redshifts. The same procedure is 
repeated for clusters at z = 0.53 and z = 0.0 (Figures ^c 
and d). Although the distributions of ei and €2 are roughly 
constant in time, individual clusters may show large changes. 
These occur both ways, i.e. ei and £2 can both increase and 
decrease. The mean change in ei between redshifts z — 1.3 
and 0.53 is 0.00, and the r-m-S- scatter around this line is 
0.16. Between redshifts z — 0.53 and the present, the mean 
change is -0.01 and the r.m-s. scatter is 0.12. For £2 these 
values are -0.01 and 0.20 between z = 1.3 and 0.53 and 0.00 
and 0.25 between z = 0.53 and z = 0.00, respectively. 

For the ilo = 0.2 scenarios, this picture does not change 
qualitatively, but quantitatively the changes are somewhat 
smaller. In Figure 11 the same plots are shown as in Fig- 
ure ^ but now for clusters in the LCDMb scenario. Low-f2o 
clusters change their ellipticities in the course of time as well, 
but the dispersions are smaller than in Figure The mean 
and r.m.s. change between z = 1.3 and 0.53 are 0.02 and 
0.08 for ei and 0.01 and 0.17 for £2, respectively. Between 
z = 0.53 and 0.0 these values are -0.03 and 0.09 for ei and 
-0.03 and 0.24 for £2, respectively. However, one should keep 
in mind that the number of clusters in the LCDMb scenario 
is (much) smaller than in the hCDM scenario, especially at 
later times. The results for the LCDMc and LCDMd sce- 
nario are similar to those for the LCDMb scenario. 
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0.25. For the LCDMb scenario, the mean and r.m.s. changes 
1 between z — 1.3 and 0.53 are -0.02 and 0.20, respectively, 
and between z = 0.53 and z = 0.0 these values are -0.03 and 
'^■^ ^ 0.17, respectively. The mean value of Tzi — Tz2, with zl > z2, 
^ is negative for all scenarios, which indicates that, on average, 
T increases with time. That is, clusters become more prolate. 
This may also be the reason why the ilo = 0.2 model clusters 
have, on average, a larger value of T. They have evolved 
further than clusters in an flo = 1.0 scenario. The r.m.s. 
scatter decreases slightly with time in most scenarios. 

Our results, that the values of ei and 62 change in time, 
while the distribution of ei is constant in time, may be ex- 
plained as follows: at a specific time new particles, that were 
not within the cluster before, enter the cluster. For an indi- 
vidual cluster, the new cluster particles may enter the cluster 
via one specific direction and change its shape. However, av- 
eraged over the cluster population, the new particles fall in 
isotropically and therefore do not change the distribution of 
ellipticities of the cluster population as a whole. 

The idea that clusters become more prolate at later 
times is consistent with the findings of Salvador-Sole & 
Solanes (1993). They concluded that the elongations of clus- 
ters are consistent with clusters being prolate and that the 
elongations are mainly produced by the tidal interactions of 
1 .0 sufficiently massive nearby clusters. 
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Figure 11. Evolution of the shape of clusters in the LCDMb 
scenario. Only clusters are used that have a mass of at least 2.22 X 
10^^/i~^Mq at both redshifts. The panels have the same meaning 
as in Figure 10. 



Combining the results for all scenarios, we conclude that 
the mean values of ei^^i — ei,z2, with zl ^ z2, are always 
consistent with zero. I.e., there is no significant change in the 
average cluster ellipticities over the time interval studied. 
The r.m.s. values range from about 0.10 (for ei and Jlo < 
1.0) to about 0.25 (for £2 and = 1-0). 

Using the number-selected cluster samples, the results 
for the changes in ei and £2 are very similar to those for the 
mass-limited samples. 

It is interesting to note that the change in e-z is always 
larger than the change in ei . This may be caused by the fact 
that £2 is always smaller than ei [Eqs. (^)], and therefore its 
determination is somewhat more difficult, as values closer 
to zero are more difficult to measure (de Theije et al. 1995). 
However, because of the rather large number of particles 
per cluster, at least lOO/fio, this cannot account for the 
whole effect. One possible explanation is the following: a 
change in b only affects £2, whereas a change in c only affects 
£1. If b and c(< b) are fixed, a change in a has a larger 
influence on £2 than on £1 . Quantitatively, if one adopts the 
mean values {£1) ~ 0.55 and (£2) ~ 0.35, one expects on 
the basis of this simple argument that A£2 = (1 — {£2))/(l — 
(£i))A£i ~ 1.4A£i. This is roughly equal to the actual values 
of A£2/A£i that are found, though the scatter in A£2/A£i 
between different scenarios is rather large (A£2/A£i « 1.3 — 
2.0). 

For the hCDM scenario the mean and r.m.s. changes in 
T between z = 1.3 and 0.53 are -0.02 and 0.29, respectively. 
Between z = 0.53 and 2: = 0.0 these values are -0.02 and 



4.5 Alignments 

Dekel et al. (1984) and West et al. (1989) concluded that the 
relative orientations of cluster major axes with the direction 
towards neighbouring clusters provide a sensitive test for the 
formation of large-scale structure in the Universe. However, 
they used 'only' about 10,000 particles in their simulations, 
so their results may be influenced by resolution effects. We 
study this question again, using our simulations in which 
we measure cluster major axes using the tensor of inertia 
method (see Section 4.4). We consider here only the align- 
ments in 3-D. The alignments that result in projected 2-D 
data will be discussed in Paper II. Three types of alignment 
of a cluster with its environment are investigated. 



4-. 5.1 Alignment with nearest neighbour 

The first type of alignment considered is that between the 
cluster major axis and the direction towards the nearest 
neighbour cluster. This is the sort of alignment that Binggeli 
(1982) found observationally in projection. It may arise, e.g., 
from the tidal force of a neighbouring cluster on the cluster 
under consideration. Likewise, the filamentary structure of 
the mass (or galaxy) distribution may induce this type of 
alignment. 

The distribution of the angular difference between the 
cluster major axis and the direction to its nearest neighbour 
is shown in Figure ^ In this Figure, only clusters with a 
mass of at least 2.22x10^* h'^ Mq are included, which have a 
nearest neighbour with a mass of at least 1.00 x lQ^'^h~^ Mq 
and which is closer than 2Qh^^ Mpc. The results are plotted 
in terms of cos(A(^) because a random distribution expresses 
itself via a uniform distribution of cos(A(jf>). 

There is some dependence of these distributions on the 
cosmological scenario. However, the differences are rather 
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Table 6. Kolmogorov-Smirnov (KS) confidence levels that the clusters in our simulations do not show an alignment efi'ect. The first two 
columns indicate the scenario and the value of trg. The quantities with the subscript '144' refer to the number-selected cluster samples, 
whereas all other quantities refer to the mass-limited samples. Columns 3 to 5 give the KS— confidence levels that the cluster major 
axis is not aligned with the direction towards its nearest neighbour. The number of clusters involved in the mass-limited samples is 
given in column 3. Column 6 to 8 give the KS— confidence levels that the cluster major axis is not aligned with the mass distribution 
within 10h~^ Mpc around it. The number of clusters involved in the mass-limited samples is given in column 6. Column 9 and 10 give 
the KS-confidence levels that the cluster major axis is not aligned with the major axis of its nearest neighbour cluster. The number of 
clusters involved in the mass— limited samples is given in column 3. See the text for more details. 



scenario 


0"8 


Ni,z 


^KS,1 




^KS, 1,144 


N2 


^KS,2 






-PkS,2,144 


^KS,3 


^KS,3,144 


SCDM 


0.46 


23 


1.00 




1.00 


78 


0.13 






0.26 X IQ-i 


1.00 


0.70 


SCDM 


0.61 


309 


0.65 




0.57 


398 


0.79 X 


10- 


6 


0.39 X 10-2 


0.11 


0.23 


LCDMb 


0.90 


12 


0.85 




0.26 X 10~i 


49 


0.54 X 


10- 


-1 


0.15 X lO-"' 


0.97 


1.00 


LCDMc 


0.90 


7 


0.85 




0.34 X 10~i 


47 


0.43 






0.58 X 10-3 


1.00 


1.00 


LCDMd 


0.90 


140 


0.32 




0.75 


179 


0.36 X 


10- 


-3 


0.92 X 10-3 


0.29 


0.20 




0.46 


4 


0.91 




1.00 


23 


0.55 






0.49 


0.98 


0.78 




0.64 


148 


0.97 




0.99 


191 


0.30 X 


10- 


-3 


0.63 X 10-2 


0.90 


1.00 


hCDM 


0.44 


4 


1.00 




1.00 


34 


0.65 X 


10- 


-1 


0.74 X 10-1 


0.70 


0.61 


hCDM 


0.65 


336 


0.26 




0.36 


425 


0.24 X 


10- 


-7 


0.23 X 10-3 


0.19 


0.82 


hCDM 


1.00 


1108 


0.66 X 10- 


-6 


0.86 X 10-2 


1168 


0.24 X 


10- 


-18 


0.11 X 10-2 


0.52 


1.00 


TCDM 


0.47 


24 


0.29 




0.66 


81 


0.19 X 


10- 


-2 


0.40 X 10-3 


1.00 


0.60 


TCDM 


0.60 


259 


0.93 X 10- 


-2 


0.85 X 10-1 


338 


0.91 X 


10- 


-7 


0.40 X 10-5 


0.20 


0.94 


HDM 


0.44 









1.00 


35 


0.63 X 


10- 


-1 


0.84 X 10-2 




0.81 


HDM 


0.66 


303 


0.26 




0.37 


458 


0.83 X 


10- 


-17 


0.33 X IQ-^ 


0.18 


0.67 


HDM 


1.00 


953 


0.20 X 10- 


-9 


0.11 X 10-3 


1063 


0.91 X 


10- 


-47 


0.44 X lO-'^ 


0.20 X 10-1 


0.37 



hard to quantify because of the small number of cluster 
pairs involved in some of the scenarios. The fourth column 
of Table |^ gives the confidence levels, derived from a KS- 
test, that the distribution of cos{/S.(f>) is consistent with a 
random distribution. The third column gives the number 
of cluster pairs involved in this analysis. All scenarios are 
consistent with a random distribution, except the hCDM, 
TCDM and HDM scenarios with the largest value of crs. This 
could be mainly due to the large number of cluster pairs in 
these scenarios. If the minimum mass limit for the parent 
cluster is decreased by a factor of two, the two S7o = 0.2 
scenarios show a significant alignment effect as well, while 
the Slo = 0.8 scenario shows a marginally significant align- 
ment effect. Apparently, the significance is mainly restricted 
by the number of cluster pairs in the sample. Changing the 
minimum mass of the nearest neighbour cluster does not in- 
fluence the above results. Relaxing the constraint that the 
nearest neighbour should be within 2Qh~^ Mpc of the par- 
ent cluster only changes the result for the HDM scenario 
with (78 = 0.66. The KS-confidence level for this scenario to 
have a random distribution of cos(A0) then becomes 0.048, 
a marginally significant alignment. 

There is a somewhat stronger cluster alignment if one 
only considers cluster pairs with distances between 10 and 
20/i-i Mpc. 

To check how the above results depend on our defini- 
tion of the cluster sample, the same analysis is done for the 
number-selected samples. The results are given in the fifth 
column of Table ^ All KS-confidence levels are very similar 
to those in column 4. The significance level is different from 
that in column 4 only for the fio = 0.2 scenarios, the hCDM 



and HDM scenarios with erg = 1.00, and the TCDM scenario 
with (jg = 0.60. For the former two scenarios, the alignment 
effect is just significant for the number-selected catalogues 
while they were not significant for the mass-limited cata- 
logues. The reason for this difference is the larger number of 
clusters in the number-selected catalogues. 

From the table it appears that the significance of the 
cluster alignment with respect to its nearest neighbour in- 
creases with (jg. In other words, the alignment effect gets 
stronger at later times in the evolution. This can be un- 
derstood by realizing that it takes some time for the tidal 
torque of a neighbouring cluster to build up this effect. 

4-5.2 Alignment with environments 

The second type of alignment we investigate is that between 
the cluster major axis and the particle distribution around 
the cluster. This is similar to the effect described by Ar- 
gyres et al. (1986) and Lambas, Groth & Peebles (1988), 
who found that galaxy counts are systematically high along 
the line defined by the projected major axis of a cluster or 
of its dominant galaxy. The effect extends to at least lbh~^ 
Mpc from the cluster centre. 

The distributions of the angular difference between the 
cluster major axis and that of the mass distribution within 
lO/j-i Mpc from the cluster are shown in Figure All 
clusters in the mass-limited cluster catalogues are used for 
this analysis. There is a significant alignment in almost all 
scenarios, showing that a cluster is strongly aligned with 
its surroundings. The alignment is usually much stronger 
than the alignment of the cluster major axis with its near- 
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cos(A0) cos(A0) 
Figure 12. Angular difference between tlie cluster major axis and 
the direction to its nearest neighbouring cluster. The distributions 
are for clusters with a mass of at least 2.22 X lO^h-lM©, while 
the nearest neighbour has a mass of at least 1.00 X 10^'^h~^ Mq 
and is within 20 Mpc. The different panels and curves have 
the same meaning as in Figure 2. 
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Figure 13. Angular difference between the cluster major axis and 
the mass distribution within 10?i~^ Mpc around it. The different 
panels and curves have the same meaning as in Figure 2. 



cance levels for the mass-limited sample are limited by the 
small number of clusters used in some cases. 



est neighbour, which generally is at a distance smaller than 
10h~^ Mpc. The seventh column in Table ^ gives the KS- 
confidence levels that the distribution of cos(Acj!>) results 
from a random orientation. Column 6 gives the number of 
clusters that is used in this analysis. For almost all f2o = 1-0 
scenarios, the significance is very high. Only the scenarios 
with the lowest values of as have lower significances. The 
LCDMc scenario does not show a significant alignment ef- 
fect, while the LCDMb scenario shows some effect. These 
results may again be somewhat misleading because these 
scenarios have fewer clusters with masses above the mass 
limit. 

If one uses all particles within 20h~^ Mpc of the parent 
cluster, the alignment signal is still present, though some- 
times somewhat less significant. If the surrounding mass dis- 
tribution in the annulus between 20 and 30h~^ Mpc around 
the cluster is used, only the HDM scenario with as — 1.00 
shows a positive detection of the alignment. 

In column 8 in Table ^ we give the results if the number- 
selected catalogues are used to do the above analysis instead 
of the mass-limited catalogues. In this case, almost all sce- 
narios show a significant alignment. Apparently, the signifi- 



4-5.3 Alignment with nearest neighbour major axis 

The third type of alignment that is evaluated is that be- 
tween the major axes of two neighbouring clusters. This 
type of alignment is, of course, highly correlated with the 
first type of alignment that we considered. However, it is 
not exactly the same. Consider two clusters which are mu- 
tual nearest neighbours. If both make an angle of, e.g., 30 
degrees with the line connecting them, their relative angle 
is between and 60 degrees. Furthermore, being a nearest 
neighbour is not a commutative property. If cluster j is the 
nearest neighbour of cluster i, but has itself cluster k ^ i 
as nearest neighbour, then this type of alignment may differ 
significantly from the first. 

As before, only clusters whose mass is larger than 2.22 x 
lO^'*/i~^Af0, with a nearest neighbour with a mass larger 
than 1.00 x 10^*/i"^Mq and closer than 20 h'^ Mpc, are 
taken into account. We find that there is no alignment of 
major axes of neighbouring clusters. The ninth column in 
Table ^ gives the KS-confidence levels for the model clusters 
to show no alignment. The number of cluster pairs involved 
in this analysis is given in column 3. Lowering the minimum 
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mass limit of the parent cluster by a factor of two does not 
change this result. 

Similarly, the results for this type of alignment do not 
change if one uses the 144 most massive groups in the simu- 
lations. So the previous results are not due to small number 
statistics. For completeness, the KS-confidence levels for the 
number-selected sample are shown in column 10 of Table 6. 



before (Section ^j), for the spectral parameter F we can 
only compare scenarios for which the value of ag differs as 
well. So those scenarios will not yield a direct indication of 
the influence of F only, but of the combination of F and 
erg. In this Section, we will only discuss the combinations of 
cosmological parameters and cluster properties that appear 
correlated. 



^.5.^ Overall alignment properties 

Summarizing the alignment properties of clusters with their 
nearest neighbouring clusters and their environments, one 
may conclude that clusters tend to be strongly aligned with 
their surrounding mass distribution in almost all scenar- 
ios. The alignment with the direction towards their nearest 
neighbour cluster is less prominent. It is only significant for 
clusters in the hCDM and HDM scenarios, which have the 
largest value of erg, and marginally significant in the LCDMb 
and LCDMc scenarios, both with fio = 0.2, and the TCDM 
scenario. The alignment between cluster major axes of near- 
est neighbours is not significant in any of the scenarios. Even 
though a cluster may be aligned with the direction towards 
its nearest neighbouring cluster, this nearest neighbour itself 
is more likely to be aligned with its immediate surroundings 
than with the former cluster. 

The different strengths of the various types of alignment 
suggest that although clusters are aligned with their nearest 
neighbour, they may be rotated around this direction with 
almost random rotation angles. Note that clusters are not 
necessarily each others nearest neighbours. For neighbours 
within 20 Mpc, there is only a 52-60% probability that 
cluster j has cluster i as its nearest neighbour if j is nearest 
neighbour of i. This probability is almost identical in all 
scenarios and causes nearest neighbour cluster major axes to 
be less aligned than the major axis of a cluster with respect 
to the direction towards its nearest neighbour. 



5 RELATIONS BETWEEN COSMOLOGICAL 
PARAMETERS AND CLUSTER 
PROPERTIES 

Our set of scenarios is chosen in such a way that one can 
compare scenarios that differ in the value of one cosmologi- 
cal parameter only. The influence of this particular param- 
eter on the cluster properties can then be investigated. One 
may even try to get scaling relations which describe the de- 
pendence of specific cluster properties on the values of the 
cosmological parameters. This analysis is complicated by the 
fact that the clusters themselves can be very different in dif- 
ferent scenarios. For example, the clusters that are identified 
in the f2o = 0.2 scenarios have a much smaller ma ss than 
clusters in the fio = 1.0 scenarios (see Section [4.l[ ). If one 
wants to investigate the influence of fio on, e.g., the cluster 
shape, one can either use the whole cluster sample or only 
the clusters in a specific mass range. In the former case, the 
correlation between cluster mass and shape (de Theije et al. 
1995) will influence the results. 

In Table ^ we summarize the parameters of which the 
effect on cluster properties are investigated. Column 2 gives 
the probed values of the particular parameter. Column 3 
shows the scenarios used for the comparison. As discussed 



5.1 R.m.s. mass fluctuation on scales of 8/i 

Mpc, (78 

When one compares all scenarios which differ only by 
the value of as, one finds that a larger value of as pro- 
duces a larger number of clusters with a mass of at least 
2.22 X lO"ft"^M0 (see Table ||), and that each cluster in- 
dividually contains more mass. These correlations are con- 
sistent with previous findings of, e.g.. White, Efstathiou & 
Frenk (1993) and Eke et al. (1996). Using the functional form 
N oc agi+'=2'"» (Eke et al. 1996), one obtains ci « 5.5 - 6.0 
for all f2o = 1-0 scenarios. The parameter C2 depends on the 
scenario and ranges from -5.1 for the HDM scenario to 1.4 
for the scenario. Note that an exponential of the form 
A*' oc exp(ca8) does not provide a reasonable fit. 

For the = 1.0 scenarios, the fraction of high— a^^g 
clusters increases with ag only for the hCDM and HDM 
scenarios, as does the largest value of aj^g. The median value 
of ajQg is constant if one restricts the analysis to clusters in 
the mass range 1.5 x lO^'^/i-^M© < M < 2.0 x lO^/i-Ufg. 
This is expected because a^^g correlates with M (see Figure 
and using a fixed and very small mass range will result 
in clusters which have almost identical values for a^^g. 

These results are independent of the mass limit that 
one applies, at least in the range 1.0 — 5.0 x 10^^/i~^Mq. 
These results cannot be checked directly for S7o < 1.0 from 
our scenarios. 

If one uses the set of the 144 most massive clus- 
ters, the median value of aj^g increases with as for 
the fio = 1.0 scenarios and is constant in time for the 
Q.0 < 1.0 scenarios. In particular, for the SCDM scenario 
one gets o^ios median °*- ''"g'^'^''^''''''^, for the hCDM sce- 
nario ajQg jng(jian o-~0 84±o.oi ^ ^-^^ HDM scenario 
''^los median (j~0 78±o.oi ^^i^ indicated errors are due to 
the fitting only. 

The particle mass of the clusters increase with ag as 
well. The fractional increase is about twice as large as for 
''^los median ' ^ ^® expected from the virial theorem. Al- 
though this is only shown for the f2o = 1.0 scenarios. White 
et al. (1993) showed that it holds for open scenarios as 
well. The median value of the cluster particle mass seems 
to increase significantly with ag only for the hCDM and 
HDM scenarios. For the hCDM scenario the dependence is 
given by M-^^^^^g^-^ oc ^-g "^^^^'^^ while for the HDM scenario 
the best-fitting relation is M.^^^-^^-^ oc ag'^^='=° '^^. The er- 
rors in the exponents are so large that the correlations are 
yet insignificant. For all other scenarios the median value 
of the cluster particle mass is nearly constant. These re- 
sults are independent of the mass limit that one applies, at 
least in the range 1.0 — 5.0 x IQ^'^h"^ Mq. If the number- 
selected cluster set is used, the results are similar to those 
for the median value of the cluster velocity dispersion. For 
the SCDM scenario M^^^:^^-^ oc al^^^''^-*^ , for the hCDM 
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Table 7. Summary of the cosmological parameters that were varied in the simulations. The first column gives the cosmological parameter. 
The second column shows the parameter values that were used for investigating the influence of the parameter. The third column gives 
the identiflcation of the scenarios that are used for the investigation. The fourth to eighth columns show the other parameters that 
describe a scenario. 
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scenario A^median Cg^^*"'^^, and for the HDM scenario 
-^median '^8 . For the low-f2o scenarios the median 

value of the cluster particle mass is constant in time. 

A larger as produces larger cluster peculiar velocities 
(see Table ^). The dispersion of the best-fitting Gaussian 
distribution to the 1-D cluster peculiar velocity distribu- 
tion scales, on average, linearly with as, with a scatter of 
about 10%. This linear relation is identical to that predicted 
by linear perturbation theory (upec oc Q^'^as; e.g., Peebles 
1993). 

The normalization as has a large effect on the align- 
ment properties of clusters. Especially the alignment of a 
cluster with its nearest neighbour and with the surrounding 
mass distribution are influenced by erg (see Table Sce- 
narios with a larger value of as contain clusters that are 
better aligned both with the direction towards their nearest 
neighbour and with their surroundings. 



5.2 Density parameter Qq 

Comparing the SCDM scenario with erg = 0.61 with the 
LCDMd scenario, it is clear that a higher value of the density 
parameter Slo will produce more clusters with a mass of at 
least 2.22 x lO^^/i'^Mg. To quantify this, we also include 
the LCDMb scenario at the epoch where as ~ 0.6. This 
is between aexp = 0.44 and aexp ~ 0.55, and the value 
of is then equal to 0.34. We will refer to this scenario as 
LCDM034 from now on. The number of clusters with a mass 
of at least 2.22 x lO^*/i~^M0 is about 23 in this scenario. 
Using the functional form oc Q'^^^'^^^ to describe the 
number of clusters with a mass of at least 2.22 x lO^^/i"^/©, 
analogous to Eke et al. (1996), one obtains ci = 1.95 and 
C2 = 2.09. An exponential of the form A'' oc exp(cr2o) does 
not provide a good fit. 

The density parameter influences the velocity dispersion 
of clusters above a certain mass. This can be seen in Figure 
[l4|a. The difference in p(> fj^g) between the SCDM and 



LCDMd scenarios is not very signiflcant, but that between 
the SCDM and LCDM034 scenarios is. The latter lacks the 
high-cr clusters that the SCDM scenario contains. If one uses 
only the 144 most massive groups in both scenarios one finds 
that the velocity dispersions of these are larger for higher Q,o, 

''"los,median °^ 

Figure ^b shows the cumulative distribution of cluster 
mass for the SCDM, LCDMd and LCDM034 scenarios for aU 
clusters in the mass-limited catalogues. All distributions are 
scaled to the same number of clusters. The SCDM scenario 
has a few clusters that have a somewhat larger mass than the 
most massive clusters in the LCDMd scenario. The curve for 
the LCDMd scenario is systematically, though only slightly, 
below that of the SCDM scenario. This is consistent with the 
findings of, e.g.. White et al. (1993). The same conclusions 
hold if one compares the SCDM scenario with the LCDM034 
scenario, but the differences are then even larger. If one uses 
the number-selected catalogues in both scenarios, the clus- 
ter mass is largest for the Q.o = 1.0 scenario (see Figure 
The median value of the cluster particle mass for clusters in 
the number-selected sample scales as M-^^^:^^-^ oc JlJ^"'*. 

In Figure ^c we show the distributions over 3-D el- 
lipticities for the SCDM, LCDMd and LCDM034 scenarios 
and for clusters above the mass limit of 2.22 x lO^/i-^M©. 
No significant differences between these distributions are de- 
tected. This may seem surprising because previous studies 
showed that the cluster ellipticity increases with f2o (Mohr 
et al. 1995, de Theije et al. 1995). However, the mass range 
involved here differ s subs tantially between the different sce- 
narios (see Section 2.2.2). As was pointed out by Struble & 
Ftaclas (1994) and de Theije et al. (1995), the more mas- 
sive clusters are less elongated than the less massive ones. 
This anti-correlation between e and M is also detected in 
the present simulations. 

To investigate whether ei and £2 are smaller for clusters 
in low-fio scenarios when clusters of the same mass are con- 
sidered, we show in Figure ^ the cumulative distributions 
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Figure 14. a: Cumulative distribution of cluster velocity dispersion for the SCDM scenario with erg = 0.61 {Qq = 1.0, solid line), the 
LCDMd scenario (SIq = 0.8, dotted line) and the LCDM034 scenario (Qq = 0.34, dashed line). Only clusters in the mass-limited samples 
are included and all three distributions have been scaled to a total cluster density of 8.6 X lQ~^h'^ Mpc~'^. b: Same as a, but for cluster 
masses, c; Distributions of cluster ellipticities. The thick lines are the distributions of ei, whereas the thin lines indicate the distributions 
of t2- Only clusters with a mass of at least 2.22 X lQ^'^h~^ Mq are used. The curves have the same meaning as in panel a. 
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narios, respectively, with dispersions of 0.10, 0.11 and 0.13. 
The mean values of 62 are 0.39, 0.38 and 0.34, respectively, 
with dispersions of 0.15, 0.16 and 0.11. Interpolating lin- 
early, we find (ei) = 0.48 + 0.13 00 and (£2) = 0.31 + 0.08^0 
for clusters in the same mass range. 

The dispersions of the fitted Gaussian distributions to 
the distributions of cluster peculiar velocity indicate that 
a somewhat smaller value of fio results in slightly smaller 
cluster peculiar velocities, consistent with previous work of, 
e.g., BahcaU et al. (1994) and Gramann et al. (1995). The 
difference in dispersions between the SCDM and LCDMd 
scenarios is 28 km s~^, or about 13%, while the difference 
between the SCDM and LCDM034 scenarios is 42 km s"^ or 
about 17%. As linear perturbation theory predicts a relation 
'^pec oc QP'^ for constant erg (e.g., Peebles 1993), we fit the 
power-law relation Cpg^ c< O'''. This does not provide a 
good fit. Neither does an exponential o-pg^, oc exp(6Slo) 
give an acceptable fit. 



Figure 15. Cumulative distributions of cluster ellipticity ei 
(thick lines) and £2 (thin lines) for the SCDM scenario (solid 
lines, 90 clusters), the LCDMd scenario (dotted lines, 227 clus- 
ters) and the LCDM034 model (dashed lines, 8 clusters). Only 
clusters that have a mass in the range (1.5 — 2.0) X lO^'*h~^M0 
are used. The mean values of ei are 0.62, 0.60 and 0.53 for the 
SCDM, LCDMd and LCDM034 scenarios, respectively 



of £1 (thick lines) and £2 (thin lines) for the SCDM scenario 
with ag = 0.61, the LCDMd scenario (having o-g = 0.60) and 
the LCDM034 scenario (having ag « 0.60). These scenarios 
differ only in the value of Q,q . Only clusters in the mass range 
1.5 X lO"/i"^M0 < M < 2.0 X 10"/i"^Mo are included. 
This is to eliminate the correlation between £i and M. It is 
clear that clusters in the same mass range are less fiattened 
in the f2o = 0.8 scenario than in the fio = 1.0 scenario, 
consistent with the previous findings of Mohr et al. (1995) 
and de Theije et al. (1995). The mean values of £1 are 0.62, 
0.60 and 0.53 for the SCDM, LCDMd and LCDM034 sce- 



5.3 Spectrum 

Comparing the SCDM with the scenario, both with 
erg = 0.46 or with erg ~ 0.61 — 0.64, one finds that the num- 
ber density of clusters is about twice as large for the SCDM 
scenario as for the scenario. The SCDM spectrum has 
more power on scales less than 8/i~^ Mpc, as its effective 
power law index in this fc-region is -1. Apparently, this ex- 
tra power on the somewhat smaller (« 2h~^ Mpc) scales 
stimulates cluster formation considerably through merging 
of smaller clumps. 

The comparison between the SCDM and scenarios 
with (Tg = 0.61—0.64 is shown in Figure p^. The line-of-sight 
velocity dispersions are slightly larger for the scenario. 
The differences in cluster mass between both scenarios are 
larger. Clusters in the k~^ scenario have larger masses than 
those in the SCDM scenario. 

The shape of the power spectrum has some influence 
on the distribution over ellipticities. The ellipticities for the 
SCDM spectrum are somewhat larger than those for the 
fe'^^-spectrum, especially for ag ~ 0.61 — 0.64. The KS- 
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Figure 16. Same as Figure The solid line is for the SCDM scenario with ag = 0.61 and the dotted line is for the k ^ scenario with 
(78 = 0.64. 



confidence levels for both distributions to be the same are 
0.19 (for £i) and 0.24 (for ea) for ag = 0.46, and 3.6 x 10"^ 
(for ei) and 1.1 x 10"^ (for £2) for erg ~ 0.61 - 0.64. The 
latter are thus significant. If one uses the 144 most massive 
groups in both scenarios, the KS-confidence levels are very 
similar to those for the mass-limited samples. 

The shape of the power spectrum has a significant effect 
on the cluster peculiar velocity. For the scenario, the 
cluster peculiar velocity is almost twice as large as for the 
SCDM scenario. This is due to the extra power of the 
scenario on the very large scales. The results for the number- 
selected samples are the same. 

The index of the primordial spectrum n only has some 
influence on the cluster number densities and peculiar ve- 
locities. For the number densities, the number of clusters 
scales as A'^ oc n'', with 5 « 1.0, although for various values 
of as the value of 5 varies from 0.8 to 1.2. The peculiar ve- 
locity is larger by about 20 km s~^ (or about 8%) for the 
TCDM scenario with n = 0.8, due to the extra large-scale 
power. Using the functional form cTpg^, 3]^ oc n~'^ one gets 
T = 0.41 ± 0.20. Because n does not correlate strongly with 
the cluster properties, its value can be chosen in order to fit 
the COBE data. Cen et al. (1992) found that n = 0.7 - 0.8 
is the most interesting range. 



5.4 Hubble parameter h 

For the SCDM and hCDM scenarios with erg ^ 0.61 - 0.65 
and for the LCDMb and LCDMc scenarios, the number of 
clusters with a mass of at least 2.22 x 10"/i"^Mo is al- 
most identical, so h does not have a large infiuence on the 
cluster number densities. Of course, there would be a large 
difference in the number of clusters if one would express 
the mass of a cluster in units of AIq instead of H'^Mq. In 
that case, the LCDMb scenario would have 201 clusters with 
M > 2.22 X W^^Mq. Qualitatively, this is expected as the 
Universe is twice as old for the LCDMb scenario as it is for 
the LCDMc scenario (see the last column of Table Us- 
ing only the comparison between the LCDMb and LCDMc 
scenario, one finds that the number of clusters with a mass 
of at least 2.22 x lO^^M© scales as cx h~°', with a in- 
creasing monotonically from 0.49 at aexp ~ 0.55 to 2.29 at 
aexp = 1.00. 



A lower value of h results in slightly larger cluster pe- 
culiar velocities (see Table |^ . This is especially true if erg is 
not too high, i.e., as < 0.6. For the dispersion of the 3-D 
peculiar velocity distribution we find the empirical relation 
(jpgg 3p oc /i"'', where /3 « 0.25 — 0.35. 

The infiuence of h on the cluster properties is similar if 
one uses the 144 most massive groups. 



5.5 Spectral parameter F 

The number of clusters with mass larger than 2.22 x 
lO^*/i~^M0 decreases with increasing F. For the LCDMc 
scenario one gets A'' oc p-i-55±0 02^ though one has to keep 
in mind that the value of erg changes as well, from 0.47 to 
0.90. The actual relation between A'^ and F may thus be 
even steeper. A similar result holds for the LCDMb scenario, 

^^p-1.43±0.01^ 

The median values of the cluster particle mass and ve- 
locity dispersion hardly change with F for the LCDMb and 
LCDMc models. 

The peculiar velocities increase with F according to 
p0.32±o.oi Jqj. ^jjg LCDMc scenario, and according to 
p0.38±o.oi fpj. ^Yie LCDMb scenario. This is surprising, as 
a larger F means that there is less large-scale power (Efs- 
tathiou et al. 1992). Furthermore, for the epochs where F 
is smaller, erg is larger. This would suggest an even steeper 
increase of the peculiar velocity towards the epochs when 
F is small. Possibly the results are infiuenced by the rather 
small number of clusters. 



6 CONCLUSIONS AND DISCUSSION 

The aim of this paper was fourfold: (1) Present the set of 
simulations which will be used in Paper II to select the 
scenario that is most consistent with the observations. (2) 
Study the intrinsic properties of clusters of galaxies for dif- 
ferent cosmological scenarios. (3) Investigate which cosmo- 
logical parameters have the largest infiuence on these cluster 
properties. (4) Obtain scaling relations between the cosmo- 
logical parameters and the cluster properties. These scaling 
relations can be used to 'interpolate' between existing sce- 
narios. 
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The following conclusions may be drawn from the 
present analysis, for the range of parameters studied and 
the scenarios studied: 

• (T8 (in the range 0.44-1.00): 

- The normalization as correlates positively with the 
cluster number density, as is expected. This is consistent 
with the earlier findings of, e.g.. White et al. (1993). 
Fitting N oc al^+'"' ""^ , we get ci « 5.5 - 6.0 for all D,o = 
1.0 scenarios, while C2 differs between various scenarios. 

- The median value of the cluster line-of-sight veloc- 
ity dispersions is almost independent of as for the r^o = 
1.0 scenarios. Only for the hCDM and HDM scenarios 
there is a slight correlation, although the errors are large. 
These results are independent of the mass limit that one 
applies, at least in the range 1.0 — 5.0 x W^'^ Mq. Us- 
ing the number-selected cluster samples, ay^^ median 
^0.78-0.90 £qj. ^jjg ^ ^ Q scenarios. 

- The median value of the cluster particle mass in- 
creases somewhat with erg for the hCDM and HDM sce- 
narios, though the errors are very large. These results are 
independent of the mass limit that one applies, at least 
in the range 1.0 - 5.0 X lO"/i"^M0. Using the number- 
selected cluster samples, M-^^^^^^^ oc a^"^^'^'^^ . 

- The value of erg does not influence the cluster ellip- 
ticities. Although for individual clusters the ellipticity 
changes, the distribution of ellipticity is independent of 
erg. 

- A larger value of erg results in significantly 
larger cluster peculiar velocities. Empirically, we found 
erpgc 2D oc erg, as is predicted by linear theory. 

- erg does affect the cluster alignment with its nearest 
neighbouring cluster and with its surroundings. The an- 
gular difference between the major axes of neighbouring 
clusters does not, or only marginally, depend on erg. 

• Qo (in the range 0.34-1.00): 

- A larger value of fio results in an increase in the 
number of clusters. Fitting the functional form A'^ oc 
Qci-Hc2no obtained Ci = 1.95 and C2 = 2.09. Chang- 
ing Qo changes the shape of the distribution of clus- 
ter line-of-sight velocity dispersions for clusters above 
a certain mass, and it results in considerably larger 
velocity dispersions for the 144 most massive clusters, 

Qa:0.39±0.08 

"los,median "o 

- A larger Slo produces more high-mass clusters both 
in absolute and in relative sense, consistent with earlier 
resuhs of White et al. (1993) and Jing & Fang (1994). 
This holds for both the mass-limited and the number- 
selected cluster catalogues. For the latter catalogues, 

^"^median °^ "o 

- A low value of fio produces more spherical clusters 

for a specific mass range. This was already concluded 
by de Theije et al. (1995) and Mohr et al. (1995). How- 
ever, for the sample of clusters above a minimum mass 
threshold or for the sample of the A'^ most massive clus- 
ters, the correlation between Qo and e disappears. This is 
because the cluster mass is very different in the different 
scenarios, and because e and M are anti-correlated. 

- A larger value of fio results in an increase in the 
cluster peculiar velocity, consistent with earlier findings 



of, e.g., Bahcall et al. (1994). However, a fit of the form 
erpgg -j^]-) oc f2g as is predicted by linear theory, where 
7 = 0.6, does not provide a good representation. 

• Spectrum (SCDM versus k~^): 

- The shape of the spectrum clearly influences the 
cluster number density. The SCDM spectrum produces 
more clusters than does the spectrum. The clusters 
in both scenarios have a similar velocity dispersion. 

- The spectrum has some influence on cluster mass 
as well. Clusters in the scenario have a larger mass 
than in the SCDM scenario. 

- The SCDM spectrum produces clusters that are 
somewhat more elongated than does the power spec- 
trum. The difference is only significant for erg > 0.6. 

- The k~^ spectrum produces cluster peculiar veloc- 
ities that are almost twice as large as for the SCDM 
spectrum. This is due to the extra power on large scales. 

- For smaller values of n, the cluster peculiar velocity 
is larger because of the extra power on large scales. 

• h (in the range 0.5-1.0): 

- The Hubble-parameter h affects the cluster peculiar 
velocity slightly: the 3-D cluster peculiar velocities scale 
as ft"'', with f3 « 0.25 — 0.35. This is because for smaller 
values of h the Universe is older. 

• F (in the range 0.20-0.43): 

- The number of clusters decreases with F according 
to oc r~-i-55±o o2 foj. the LCDMc scenario, though 
one should keep in mind that the different scenarios that 
are used do also have a different value of erg. 

- The cluster peculiar velocity scales as p~0 32±o.oi 
for the LCDMc scenario and as r-"'^***"^"^ for the 
LCDMb scenario. This is quite surprising as a larger 
F indicates less large-scale power. 

In summary, we conclude that erg has the largest influ- 
ence on the cluster properties. This is not surprising because 
the mass within the virial radius of a rich cluster is very close 
to the mass enclosed within a sphere of 8h~^ Mpc in the un- 
perturbed Universe (e.g., Evrard 1989). Almost all cluster 
properties change if erg is varied, fio has a large impact on 
the cluster number density, mass, and peculiar velocity. In 
addition, for the number-selected cluster set, relatively more 
high-erjgg clusters are expected if Qo is larger. More power 
on larger scales produces somewhat more elongated clusters 
and larger cluster peculiar velocities. The other parameters, 
the spectrum and h, correlate less strongly with the cluster 
properties. 

The cluster peculiar velocity is the cluster property that 
depends on the largest number of parameters of the fluctu- 
ation scenario. Two difficulties affect the determination of 
the most consistent cosmological scenario purely on the ba- 
sis of cluster peculiar velocities. First, the large sensitivity 
of cluster peculiar velocities on all cosmological parameters 
make it very hard to disentangle these parameters and de- 
termine each of them separately. Fortunately, quite a few 
relations between cosmological parameters are known from 
other studies. E.g., Big Bang Nucleosynthesis and the fact 
that clusters cannot consist of more than 100% baryons put 
severe limits on the Hubble parameter /i as a function of Qo 
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(David, Jones & Forman 1995). Secondly, cluster peculiar 
velocities are very hard to determine observationally. Mould 
et at (1991, 1993), e.g., quoted errors of 300 to 800 km 
s~^. Different studies sometimes yielded very different val- 
ues for vpec for the same cluster. Very recently, Giovanelli 
et al. (1997) obtained the peculiar velocities for a sample 
of 22 groups and clusters. Although the sample is rather 
small, the uncertainties are considerably smaller (about 150 
km s~^) than those in previous studies and it will be very 
worthwhile to extend this data set to a larger number of 
clusters. 

As the cluster number density depends on the normal- 
ization (Tg, on the density parameter Qq and on the shape of 
the power spectrum, it can be used to discriminate between 
different scenarios (White et al. 1993, Eke et al. 1996). How- 
ever, the cluster number density in the simulations depends 
on the exact definition of a cluster. The most straightfor- 
ward definition is, of course, to apply a mass threshold and 
consider all objects with a mass larger than this threshold to 
be clusters. But it is hard to get a reliable mass estimate of 
a cluster from galaxy positions and velocities. Better mass 
estimates may be obtained from X-ray measurements and 
gravitational lensing. Another way to define a cluster may 
be to use the cluster line-of-sight velocity dispersion and use 
only those objects which have a aj^g larger than some well- 
chosen value. With the large redshift surveys coming up in 
the very near future, and with the ENACS-survey (Katgert 
et al. 1996) already being completed, this may prove to be a 
suitable manner to define a cluster. However, one does not 
necessarily pick out the most massive objects because of the 
scatter in the crj^g — M-relation (see Figure ^). 

It is promising to note that most of the cluster prop- 
erties discussed do not depend critically on the cluster def- 
inition. That is, most results are qualitatively the same for 
the mass-limited and the number-selected cluster sample. 
However, when comparing to observations one should try to 
construct a model cluster catalogue that is complete in, e.g., 
richness or X-ray temperature. 

In Paper II, we will compare the properties of the 144 
most massive clusters in the simulations with many obser- 
vations. We then try mimic as closely as possible the way 
in which real clusters are observed. The results of the com- 
parison with observations, together with the results of other 
studies, can then be used to find the cosmological scenario 
which is most consistent with all cluster properties. If nec- 
essary, the scaling relations can be used to 'interpolate' be- 
tween scenarios. For the best-fitting scenario, a higher res- 
olution cluster catalogue can then be constructed for which 
each cluster will be simulated individually at high resolution. 
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